Electronic structure and light-induced conductivity in a transparent refractory oxide 
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Combined first-principles and experimental investigations reveal the underlying mechanism re- 
sponsible for a drastic change of the conductivity (by 10 orders of magnitude) following hydrogen 
annealing and UV-irradiation in a transparent oxide, 12CaO-7Al2C>3, found by Hayashi et al. The 
charge transport associated with photo-excitation of an electron from H, occurs by electron hopping. 
We identify the atoms participating in the hops, determine the exact paths for the carrier migra- 
tion, estimate the temperature behavior of the hopping transport and predict a way to enhance the 
conductivity by specific doping. 
PACS numbers: 71.20.-b, 72.20.Ee, 72.40. +w, 81.05. Je 



The rare combination of transparency and conductiv- 
ity has a wide range of technological applications 1 . Until 
recently, optimum optical transmission and useful elec- 
trical conductivity (>10 3 S cm -1 ) was attained by the 
introduction of selected dopants into a wide-bandgap ox- 
ide. Hence the excitement surrounding the report that 
Hayashi et al found a new way to convert a transparent 
oxide into a persistent electronic conductor 2 with the po- 
tential for less expensive and more flexible optoelectronic 
devices (including flat panel displays, solar cells, invisible 
circuits) as well as three-dimensional holographic mem- 
ories. Despite careful experimental studies, no definitive 
understanding has been reached on the underlying mech- 
anism responsible for this new dramatic effect. In this 
combined theoretical and experimental investigation, we 
address the nature of this novel phenomenon which is 
crucial for further progress and make predictions on how 
to improve the conductivity toward the level where it 
becomes useful for practical applications. 

We have reproduced in bulk form the single crystal 
and thin film results of Hayashi et al by solid state re- 
action of high purity oxides 3 (99.99+%), followed by hy- 
drogen treatment (5%H 2 /95%N 2 ) of sintered pellets at 
1,300 °C for 2 h, and by quenching to room temperature. 
UV-treatment was performed under a mercury arc lamp 
(40 min at 20 mW cm~ 2 and a wavelength of 220 nm). 
This turned a surface layer of *~16 [im thickness from 
white to green, with a conductivity changed from 10~ 10 
to 0.6 S cm -1 (by four-point method). We also irradi- 
ated loose H-treated powders, which also turned from 
white to a uniform green color, and measured their con- 
ductivity by a "powder-solution-composite" method 4 to 
be 0.41 S cm -1 . Both values agree with those in Ref. 2 
and provide confidence in our experimental confirmation 
of the predictions reported below. 

From the theoretical side, we demonstrate with den- 
sity functional calculations that the conductivity appears 
to result from electrons excited off the hydrogen ions 
into conduction bands formed from Ca 3d states. The 
changes found in the electronic properties caused by the 
UV-irradiation of H-doped mayenite allow us to conclude 
that the observed absorption peak at 0.4 eV has a fun- 
damentally different nature than the F + center energy 



level originally assumed in Ref. 2 and followed in recent 
theoretical work 5 . Our calculated charge density shows 
clearly that electron trapping on a vacancy (i.e., forma- 
tion of the F + center) does not occur in the case of UV- 
irradiated H-doped mayenite. Instead, we identify the 
atoms participating in hopping transport; estimates of 
the characteristic parameters that determine the behav- 
ior of the hopping conductivity at low (Efros-Shklovskii 
regime) and high (Mott) temperatures agree reasonably 
well with experiment. Finally, we found a strong de- 
pendence of the electronic transport on the particular 
hopping centers and their spatial arrangement that al- 
lows the prediction of rates of temperature treatment and 
dopants which will enhance the conductivity. 
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FIG. 1. (a) Three of the 12 cages constituting a unit cell 
of mayenite. An O 2- anion (abbreviated as O*) is located 
inside one of six cages so that the distances between O* and 
its six neighbor Ca atoms vary from 2.36 to 3.95 A. (b) Total 
and partial DOS of 12CaO-7Al 2 3 ; E F is at eV. The Al 3d 
states form the band centered at ~28 eV. 



The crystal structure of 12CaO-7Al203 has the space 
group /43c?, with two formula units in the unit cell and 
a lattice parameter of 11.98 A (Ref. 6, 7). This frame- 
work includes 64 of the oxygen atoms; the remaining two 
oxygen ions (abbreviated as O* hereafter) which provide 
charge neutrality 7 are located inside the cages, Fig. 1(a). 
These two O* are distributed between 24 sites that pro- 
duce a structural disorder. From full-potential linearized 
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augmented plane-wave 8 (FLAPW) total energy calcula- 
tions of seven structures with different distances between 
these two oxygen ions (ranging from 5.1 A to 10.4 A), 
we found that the 0*s tend to maximize the distance 
between them by forming a bcc lattice with the crystal- 
lographic basis oriented randomly with respect to that 
of the whole crystal. After hydrogen annealing, changes 
in the crystal structure are not observed 9 . The incorpo- 
ration of H ions corresponds to the following chemical 
reaction: O 2- (cage) + H 2 (atm.) — > OH~ (cage) + 
H~ (in another cage), so the unit cell contains two cages 
occupied by 0*H, another two by hydrogen ions (abbre- 
viated as H* hereafter) and the remaining 8 cages are 
empty 10 . From the electronic structure calculations, we 
found that the structure with two H - ions in cages and 
released O2 gas (cf., Ref. 2) leads to an insulating state 
after photo-activation. 
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2. (a) Total DOS of H-bearing mayenite. The band 
g , is between O 2p and Ca 3d states. The DOS of 
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the impurity bands arising from the OH - complex and H* 
is enlarged by a factor of five. The photo-activation of the 
H* corresponds to the transition from H* Is states to the 3d 
states of its closest Ca (dashed line), with energy E e . (b) 
Total DOS of UV-irradiated H-doped mayenite and the joint 
DOS (thick line) calculated using the optical selection rules 
and only nearest neighbor transitions which have higher prob- 
ability, (c) Enlarged total DOS of UV-irradiated H-doped 
mayenite and a theoretical fit to the DOS near Ef (dashed 
line), see text, (d) The optimal path for electron migration 
throughout the crystal. The cube represents a full unit cell 
consisting of 12 cages (118 atoms); only atoms which give the 
dominant contribution to the band crossing Ef are shown - 
Ca (largest spheres), O* (between Ca(2) and Ca(3)) and H* 
(smallest spheres). 

The calculated 11 electronic density of states (DOS) be- 
fore H-doping is shown in Fig. 1(b): oxygen 2p states 
form the top of the valence band between —6 eV and —2 
eV (the Fermi level is taken as zero), while the bottom 
of the conduction band is formed from Ca 3d states. Lo- 
cated inside one of the six cages, O* gives a fully occupied 
impurity peak below Ep. The calculated band gap, 4.8 



eV, underestimates the experimental absorption edge 2 by 
4 %. The incorporation of H ions into the cages of insulat- 
ing mayenite results in the appearance of new bands (Fig. 
2(a)): (i) since H shares its electron with O* (and forms 
an OH - complex), filled a and non-bonding 7r bands and 
unoccupied a' band are introduced; and (ii) the Is states 
of H* form the fully occupied impurity band below Ep-- 
After the H-doping, the energy of the transitions from O 
2p states to Ca 3d states is decreased by 0.7 eV (from 
5.9 eV to 5.2 eV), in agreement with experiment 2 . The 
maximum efficiency of the photo-activation of the H*, 
which corresponds to transitions from H* Is states to 
the 3d states of its closest Ca neighbors, occurs at about 
the same photon energy (Fig. 2(a)). These results are 
in contrast with the conclusion of Ref. 2 that hydrogen 
controls the absorption edge. 

For UV-irradiated H-bearing mayenite, the electron 
excitation from H* is calculated with the model system 
where the electron is transferred equally to the d states 
of six Ca neighbors 12 . The self-consistent DOS shows the 
formation of a new band (Fig. 2(b)) with the same total 
number of states, m, as in the band gap of the system be- 
fore UV irradiation. However, since m is now two times 
larger than the number of available electrons, Ep passes 
through the middle of the band and the system becomes 
conducting. The calculated positions of the characteris- 
tic absorption bands (0.38 eV and 3.5 eV, Fig. 2(b)) are 
in good agreement with experiment (0.4 eV and 2.8 eV), 
as is the intense charge transfer transition from oxygen 
2p states to Ca 3d states which occurs at ~4.8 eV. 
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FIG. 3. Contour map of the charge density distribution 
within a slice passing through the center of (a) a cage with 
H* and (b) an empty one (vacancy). 

As shown in Fig. 2(c), the non-zero DOS at E^, g(E F ), 
is determined mainly by H*, O* and Ca states - atoms 
that are spatially well separated from each other (Fig. 
2(d)), which points to the hopping nature of the induced 
conductivity. It is striking that the arrangement of the 
atoms corresponds to the shortest electron hopping path: 
only the closest Ca neighbors of O*, Ca(2) and Ca(3) (at 
distances of 3.3 and 2.4 A, respectively), and the closest 
Ca neighbors of H*, Ca(l) and Ca(4) (both at 2.8 A), 
give significant contributions to g(Ep) (Table I). The 
longest segment on the hopping path corresponds to the 
distance between Ca(l), Ca(2) and O*, which explains 
the higher probability to find the electron on Ca(l) and 
Ca(2). Finally, the charge density distribution calculated 
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in the energy window of 25 meV below Ep, clearly shows 
the connected electron density maxima along the hopping 
path (cf., Fig. 3(a)). The charge density in the empty 
cages is insignificant (cf., Fig. 3(b)), demonstrating that 
trapping of the electron on the vacancy does not occur 13 . 

Hopping transport in systems with localized electronic 
states has been studied analytically, numerically and 
experimentally 14 " 16 . At high temperatures (T) the hop- 
ping occurs between nearest neighbors and has an acti- 
vated nature. At lower T, the activation energy for hops 
between localized centers can be reduced by enlarging the 
hopping distance - which leads to the variable-range hop- 
ping described by Mott 14 . At yet lower T, the electron- 
electron repulsive interaction between charge carriers 
may result in the formation of a Coulomb gap centered 
at Ep. The DOS at Ej? may not turn exactly to zero, 
corresponding to a 'soft' Coulomb gap 15 : this residual 
DOS at Ef is caused by correlated simultaneous hops of 
several electrons (multi-electron hopping) which plays a 
dominant role 17 at extremely low T. 

Similar to the above qualitative picture, we predict 
the behavior of the conductivity at finite T based on the 
calculated DOS in the vicinity of Ep for the irradiated 
H-doped mayenite. Clearly, this exhibits in Fig. 2(c) 
a soft Coulomb gap near Ep and, similar to Fig. 10.1 
in Ref. 15, has two peaks. These results conform to the 
Efros-Shklovskii (ES) behavior 15 of the conductivity a at 
low T: 



cr = a exp[-(T ES /T) 1 / 2 }, 



T E s = (3Ese 2 /k B k£, 



(1) 



where £ is the localization length, e - the elementary 
charge, ks - the Boltzmann constant, and k=eoe r is the 
dielectric constant (sq - the permittivity constant and 
e r - the relative dielectric constant). When T increases, 
the Coulomb gap fills. According to Ref. 18, a substantial 
DOS should be observed already at T~0.05Ae, where Ae 
is the gap width. At this T, the changeover to the Mott 
behavior 14 should occur: 

a = a exp[-(T M /T) 1 /% T M = M /k B g Q ? . (2) 

The double peak structure completely disappears at T~ 
0.4Ae. Using Ae~0.3 eV (Fig. 2(c)), we obtain the 
crossover temperature between Mott and ES behavior 
to be T c ~170 K. Indeed, for H-doped UV-irradiated 
mayenite, the J 1-1 / 4 dependence of log(a) was found 
experimentally 2 for 50<T<300 K. Further, the calcu- 
lated DOS near E F fits well (Fig. 2(c)) with the the- 
oretical one of the form g(E)—aE 2 E 2 /(E 2 + E 2 ), where 
a=(3/-7r)(A; 3 /e 6 ). In the high energy limit, the DOS has 
the constant value of g =aE 2 that corresponds to Mott 
behavior. In the opposite limit, it approaches a parabolic 
dependence (ES) which has only one physical parameter 
- the refractive index, n. We found n~2.1, while the 
experiment 20 gives n=1.6. From the fit, we obtain go— 8 
states/eV-cell and a~800 which allows us to estimate 



the characteristic T: using /3 M =7.6 (Ref. 19), (3 E s=2.8 
(Ref. 15) and the localization length £=1.5 A (which is 
comparable to the ionic radii of the atoms involved in the 
hops), we obtain Tm— 5.6x 10 6 K which agrees fairly well 
with the experimental value 2 , Tm=2x10 6 K, and predict 
that T£s~2.7xl0 4 K. Thus, our band structure calcula- 
tions give reasonable order of magnitude estimates for the 
characteristic T of the hopping transport in mayenite. Of 
course, more precise determinations require temperature 
and time-dependent calculations. 
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FIG. 4. Measured ct-T characteristics of the UV-irradiated 
H-doped 12CaO-7Al 2 3 bulk samples of 46% density. Due 
to the realization of the optimal arrangement of the hopping 
centers during the heat treatment (A— >B), the conductivity 
is enhanced after cooling the sample (B— >C). Subsequent T 
cycling (B<-»C) does not result in any changes of cr, i.e., the 
system is reversible. Cooling from any point above the con- 
ductivity maximum, as in path (D— >E), shows reduced a un- 
changed by UV irradiation at room T (E<-»D), as expected 
from the release of H. 

Since H* may also occupy cages located farther away 
from the one with an 0*H, we investigated this arrange- 
ment and found it to be less favorable energetically - by 
96 meV. However, under the rapid cooling of the sample 
annealed at 1,300 °C some H* can become "frozen" into 
these positions. After UV irradiation, the corresponding 
hopping path involves a very long segment between two 
Ca atoms (now increased from 3.7 to 6.5 A), so that the 
DOS at Ep becomes zero and the hopping has an acti- 
vated behavior. If now T is increased to a level sufficient 
for hydrogen migration, these H* will be able to diffuse 
into the energetically favorable positions which facilitates 
the hopping and leads to an increase of the conductivity. 
We confirmed these predictions experimentally, as shown 
in Fig. 4. Moreover, our samples showed non-reversible 
behavior once T is increased above T m , the conductivity 
maximum: the decay of the conductivity is associated 
with hydrogen release from the sample which breaks up 
the hopping path. 

If one focusses on the strong dependence of the photo- 
induced conductivity on the particular atoms participat- 
ing in the hopping and their spatial arrangement, one is 
able to predict the possibility of varying the conductivity 
by proper doping. For example, we expect that Mg sub- 
stitution can lead to a decrease in the conductivity once 
Mg substitutes one of the Ca atoms involved in the hop- 
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ping, since its 3d states will lie much higher in energy 
than that of Ca, while Sr substitution may not result 
in any significant changes in the conductivity. Instead, 
one way to possibly improve the electronic transport of 
the material consists in increasing the concentration of 
hopping centers. For this, the poorly bonded 2 ~ ion 
inside a cage of mayenite, is easily replaced by 2F~ or 
2C1~ ions that will allow some enhancement of the con- 
ductivity since the number of hopping centers and the 
number of carriers is now doubled. Recent observations 
agree with our model and its predictions: proton implan- 
tation was found to increas the conductivity as compared 
with thermal H-treatment 21 . Our additional calculations 
show that along with the increased concentration of the 
encaged electrons, the proton implantation will result in 
the appearance of new unoccupied states in the band 
gap making H + one of the hopping centers that creates 
a conductivity channel and so enhances the transport. 

Finally, the encouraging findings obtained in this first 
ab-initio study of hopping transport suggest that impor- 
tant results may also be obtained when a similar ap- 
proach is applied to such systems as doped semiconduc- 
tors, manganite perovskites, carbon nanotubes, etc 22 . 
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TABLE I. The relative contributions per atom to the DOS 
at Ep. Similar percentages are obtained from the integrated 
DOS in the energy range of E F ±0.025 eV. 
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